Search for: All records

We conduct three-dimensional direct numerical simulations to investigate the mixing, entrainment and energy budgets of gravity currents emerging from two-layer stratified locks. Depending on the density and layer thickness ratios, we find that either the upper layer or lower layer fluid can propagate faster, and that the density structure of the overall gravity current can range from strongly stratified to near-complete mixing. We furthermore observe that intermediate values of the density ratio can maximise mixing between the gravity current layers. Based on the vorticity budget, we propose a theoretical model for predicting the overall gravity current height, along with the front velocity of the two layers, for situations in which the lower layer moves faster than the upper layer. The model identifies the role of the height and thickness ratios in determining the velocity structure of the current, and it clarifies the dynamics of the ambient counter-current. A detailed analysis of the energy budget quantifies the conversion of potential into kinetic energy as a function of the governing parameters, along with the energy transfer between the different layers of the gravity current and the ambient fluid. Depending on the values of the density and layer thickness ratios, we find that the lower lock layer can gain or lose energy, whereas the upper layer always loses energy.

 

Debris flows are dense and fast-moving complex suspensions of soil and water that threaten lives and infrastructure. Assessing the hazard potential of debris flows requires predicting yield and flow behavior. Reported measurements of rheology for debris flow slurries are highly variable and sometimes contradictory due to heterogeneity in particle composition and volume fraction ( ϕ ) and also inconsistent measurement methods. Here we examine the composition and flow behavior of source materials that formed the postwildfire debris flows in Montecito, CA, in 2018, for a wide range of ϕ that encapsulates debris flow formation by overland flow. We find that shear viscosity and yield stress are controlled by the distance from jamming, Δ ϕ = ϕ m − ϕ , where the jamming fraction ϕ m is a material parameter that depends on grain size polydispersity and friction. By rescaling shear and viscous stresses to account for these effects, the data collapse onto a simple nondimensional flow curve indicative of a Bingham plastic (viscoplastic) fluid. Given the highly nonlinear dependence of rheology on Δ ϕ , our findings suggest that determining the jamming fraction for natural materials will significantly improve flow models for geophysical suspensions such as hyperconcentrated flows and debris flows. 

We investigate the submerged collapse of weakly polydisperse, loosely packed cohesive granular columns, as a function of aspect ratio and cohesive force strength, via grain-resolving direct numerical simulations. The cohesive forces act to prevent the detachment of individual particles from the main body of the collapsing column, reduce its front velocity, and yield a shorter and thicker final deposit. All of these effects can be captured accurately across a broad range of parameters by piecewise power-law relationships. The cohesive forces reduce significantly the amount of available potential energy released by the particles. For shallow columns, the particle and fluid kinetic energy decreases for stronger cohesion. For tall columns, on the other hand, moderate cohesive forces increase the maximum particle kinetic energy, since they accelerate the initial free-fall of the upper column section. Only for larger cohesive forces does the peak kinetic energy of the particles decrease. Computational particle tracking indicates that the cohesive forces reduce the mixing of particles within the collapsing column, and it identifies the regions of origin of those particles that travel the farthest. The simulations demonstrate that cohesion promotes aggregation and the formation of aggregates. Furthermore, they provide complete information on the temporally and spatially evolving network of cohesive and direct contact force bonds. While the normal contact forces are aligned primarily in the vertical direction, the cohesive bonds adjust their preferred spatial orientation throughout the collapse. They result in a net macroscopic stress that counteracts deformation and slows the spreading of the advancing particle front.

 

Emerging technologies such as deep-sea mining and geoengineering pose fundamentally new questions regarding the dynamics of gravity currents. Such activities can continuously release dense sediment plumes from moving locations, thereafter propagating as gravity currents. Here, we present the results of idealized numerical simulations of this novel configuration, and investigate the propagation of a gravity current that results from a moving source of buoyancy, as a function of the ratio of source speed to buoyancy velocity. We show that above a certain value of this ratio, the flow enters a supercritical regime in which the source moves more rapidly than the generated current, resulting in a statistically steady state in the reference frame of the moving source. Once in the supercritical regime, the current goes through a second transition beyond which fluid in the head of the current moves approximately in the direction normal to the direction of motion of the source, and the time evolution of the front in the lateral direction is well described by an equivalent constant volume lock-release gravity current. We use our findings to gain insight into the propagation of sediment plumes released by deep-sea mining collector vehicles, and present proof-of-concept tow-tank laboratory experiments of a model deep-sea mining collector discharging dense dyed fluid in its wake. The experiments reveal the formation a wedge-shaped gravity current front which narrows as the ratio of collector-to-buoyancy velocity increases. The time-averaged front position shows good agreement with the results of the numerical model in the supercritical regime. 

We investigate the interaction of a downslope gravity current with an internal wave propagating along a two-layer density jump. Direct numerical simulations confirm earlier experimental findings of a reduced gravity current mass flux, as well as the partial removal of the gravity current head from its body by large-amplitude waves (Hogg et al. , Environ. Fluid Mech. , vol. 18 (2), 2018, pp. 383–394). The current is observed to split into an intrusion of diluted fluid that propagates along the interface and a hyperpycnal current that continues to move downslope. The simulations provide detailed quantitative information on the energy budget components and the mixing dynamics of the current–wave interaction, which demonstrates the existence of two distinct parameter regimes. Small-amplitude waves affect the current in a largely transient fashion, so that the post-interaction properties of the current approach those in the absence of a wave. Large-amplitude waves, on the other hand, perform a sufficiently large amount of work on the gravity current fluid so as to modify its properties over the long term. The ‘decapitation’ of the current by large waves, along with the associated formation of an upslope current, enhance both viscous dissipation and irreversible mixing, thereby strongly reducing the available potential energy of the flow. 

We employ direct numerical simulations in order to analyze the role of double‐diffusive salt fingering in halite precipitation from hypersaline lakes. Guided by field observations from the Dead Sea, which represents the only modern deep stratified lake that precipitates halite under hydrological crisis, we consider a saturated layer of warm, salty brine (epilimnion) overlying a layer of colder, less salty brine (hypolimnion) that is also saturated. The double‐diffusive instability originating in the metalimnion gives rise to an asymmetrical pattern of less salty ascending fingers, accompanied by descending salt fingers that lose heat as they propagate through the metalimnion. The net result is a strong, downward salinity flux that leads to the undersaturation of the epilimnion, while the hypolimnion becomes oversaturated and precipitates halite. These interfacial processes within deep, hypersaline water columns in warm and dry regions suggest a potential route toward the formation of thick halite layers found in the geological record.